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Abstract. We consider the generic problem of suddenly changing the geometry 
of an integrable, one-dimensional many-body quantum system. We show how the 
physics of an initial quantum state released into a bigger system can be completely 
described within the framework of the Algebraic Bethe Ansatz, by providing an exact 
decomposition of the initial state into the eigenstate basis of the system after such 
a geometric quench. Our results, applicable to a large class of models including 
the Lieb-Liniger gas and Heisenberg spin chains, thus offer a reliable framework 
for the calculation of time-dependent expectation values and correlations in this 
nonequilibrium situation. 
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Integrable models have provided a wealth of valuable information on the equilibrium 
thermodynamics of strongly correlated systems during the last decades (see [H |2], [3] and 
references therein). More recently, much progress has also been made on their dynamical 
properties. On the other hand, the field of out-of-equilibrium dynamics is only now 
starting to become accessible to comparably exact methods [U El El [3 E] , in particular 
in the study of so-called quantum quenches (see [9]- [23] for example). In those studies 
the quench is defined as a sudden change of a local parameter of the Hamiltonian, 
typically the interaction strength or an external field. In this short note, we consider 
the case where a global parameter is abruptly changed, namely the size of the system 
before and after the quench. We baptize this the geometric quench. Earlier studies of 
such a setup, using numerical techniques, include 125]. 

In this paper, we highlight an interesting feature of the wavefunctions of integrable 
one-dimensional models of quantum many-body physics, opening the door to a large 
class of possible calculations of nonequilibrium dynamical effects in these systems upon a 
sudden change of their geometry. We start by a simple explanation using the coordinate 
representation of wavefunctions, and provide an algebraic version of the reasoning 
thereafter. The paper ends with some perspectives and suggestions for further work. 

1. Coordinate Bethe Ansatz 

We consider a quantum system initially defined within a finite space interval < x < Li 
(for clarity, we use notations appropriate for a model defined on the continuum, such 
as the Lieb-Liniger model [26] ; our reasoning also applies to lattice models such as spin 
chains by trivial modifications). We suppose that the initial Hamiltonian can be written 
as the integral of a Hamiltonian density 

Hi= I ' dxn{x) (1) 

Jo 

and that specific boundary conditions {e.g. periodic or twisted; the reasoning for the 
open case is slightly different, and is discussed in the next section) labeled by ai are 
imposed on the wavefunctions. In the case of integrable systems, a complete set of 
orthonormalizable eigenfunctions is explicitly given in each A^-particle sector by a Bethe 
Ansatz; individual wavefunctions 

^^'\{m>^}al), 0<x,<L, (2) 

are labeled by their cardinality-A set {A}^J of rapidities, solution to the Bethe equations 
associated to quantization in the space < x < Li under boundary conditions 
«!. Completeness ensures that the number of allowable solutions coincides with the 
dimensionality of the Hilbert space. One important observation, reflected in our 
notation, is that the eigenfunctions ^^dxjKA}) do not explicitly depend on the length 
Li or boundary conditions ai; rather, these dependencies are carried only implicitly by 
the sets of rapidities. This simple observation forms the basis of what we propose here. 
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The generic experiment we have in mind can be described as follows. Suppose that 
we initially prepare the system in a specific eigenstate t=0~) = ^'•^Hl^lK-^laj)- 

At time t = 0, we glue a new, empty region onto the original system, thereby changing 
its size from Li to L2 > Li, a process we call a geometric quench. From this point 
onwards, to control the subsequent time evolution exactly, all calculations must be done 
in the basis of eigenstates of the new Hamiltonian 

H2= ! ' dxn{x) (3) 







defined by integrating the same Hamiltonian density as before (so all interaction 
parameters remain unchanged) over the larger space interval up to L2. The eigenstates 
of H2 (choosing boundary conditions labeled by 0:2) are functions of the form 

^^'HWIMS)' 0<x.<L2, (4) 

where the sets of rapidities {/ij^^ ^'^^ solve the Bethe equations associated to 
quantization in the space < x < L2 under boundary conditions a2- 

At first sight, the gluing of the wavefunctions before and after the geometric quench 
is not simple. This gluing condition can be written as the nonlinear condition 

^^^'^^''"^ ^" \ if3x,G[Li,L2[. 

Decomposing this into the basis of eigenstates of H2 can now be achieved by relying on 
the crucial observation made above, namely that the wavefunctions of Hi (respectively 
H2) do not depend explicitly on Li, ai (respectively L2, 02) • This allows us to make the 
identifications 

^^'nWI{/^}£)|o<..<L, = ^«({4IWS)lo<..<L,. (6) 

In other words, in the region of space where the initial wavefunction is nonvanishing, we 
can exactly match the functional form of Bethe wavefunctions of H2 when restricted to 
the evaluation domain < Xi < Li, to the wavefunctions of Hi but this time evaluated 
using the rapidities of the state. The overlap of the initial state with a given 
eigenfunction of H2 is thus given by the well-known Slavnov formula [28] (associated to 
the original domain). 



\/(«<"({A}S)l*<'>({A}il)>('i'P'({f'}S)l*'''({A}S)) 

The exact decomposition of our initial state onto the wavefunctions of H2 is simply 
written as 

n{x},t = o^) = J2 WS)^^'^(WIMS)- (8) 

{m}^2 

The subsequent time evolution is now trivial to build back in, since are exact 
eigenstates of H2. This simple expression is the central result of our paper. 
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2. Algebraic Bethe Ansatz 

2.1. Geometric quench overlap: Periodic boundary conditions 

From the previous analysis we clearly see that the the dynamics of the geometric quench 
problem is encoded in the overlaps Fi defined by The natural framework for 

computing norms and overlaps of Bethe eigenstates is the Algebraic Bethe Ansatz. This 
approach is based on the construction of a monodromy matrix T(A), which leads to 
an algebraic formulation of the problem, see for instance [2]. The monodromy matrix, 
T G Vo ® H with Vo an auxiliary space and % the Hilbert space, should satisfy the 
intertwining relation 

R{X - fi) (T(A) ® T(^)) = (T(^) ® T(A)) R{X - /i) (9) 

where R{X) G Vq ® Vq is the i?— matrix satisfying the Yang-Baxter equation. One 
can define the transfer matrix as T(A) = TroT(A), where the trace is taken over the 
auxiliary space. The transfer matrix now acts as a generating function for a set of 
mutually commuting conserved charges, including the Hamiltonian. For simplicity we 
focus in this discussion on XXZ-type models (which includes the Lieb-Liniger gas) given 
by an i?-matrix of the form 

f /(A) \ 

g{X) h{X) 

h{X) giX) 

V /(A) / 

For this type of i?- matrix the monodromy matrix T{X) is a 2 x 2 matrix in auxiliary 
space, 

Tt\) - ( -^'^^ ^ rill 

The commutation relations of the operators A{X), B{X), C(A), D{X) G "H follow from (j9]). 
Using these operators, eigenstates of the transfer matrix (and thus of the Hamiltonian) 
can be constructed as 

N N 

(v[/|=(oinc(A.) i^)=n^(^^)io) (12) 

i=i i=i 
were the rapidities {Xj} should satisfy the Bethe equations. The pseudovacuum |0) 
((0| = lO)"'') has the following properties: 

(0|A(A) = a(A)(0| A(A)|0) =a(A)|0) (13) 
{0\D{X) = d{X){0\ D{X)\0) = d{X)\0) (14) 
(0|5(A) = C(A)|0) = 0. (15) 

The form of the vacuum eigenvalues a(A) and d{X) depend on the specific representation 
of the i?-matrix. 



RiX) 



(10) 
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In order to describe the geometric quench we decompose the Hilbert space into two 
spatial intervals: one of length Li and one of length L2 — Li (the part that is added): 
l-i = 1-Li ® 1-L2- The operators before the quench act on "Hi, while the operators after 
the quench act on "H. To be able to make a connection between the two we consider the 
generalized two-site model introduced in [27], namely the monodromy matrix T(A) of 
any Algebraic Bethe Ansatz solvable model can be written as a matrix product of two 
factors 

T(A) = r2(A)Ti(A) (16) 

with Tj(A) G Vq ® ^j- From the decomposition of the monodromoy matrix follows that 

5(A) = A2(A)fii(A) + 52(A)Di(A) (17) 

and similar expressions for the A(A), C(A) and -D(A) operators. We can write the initial 
A^-particle state (at t = 0"*^) as 

(^^^H{A})|= (^i(0|n^i(A,)j® 2(01, (18) 

while eigenstates in "H can be expressed as 

= ^Jj5i(^j)|0)ij ® [x\A2{iJij)\Q)^ + terms without overlap. (19) 

The terms without overlap are terms where not all -Bi(/i) act on |0)i and will therefore 
not contribute to the decomposition of the initial state. The overlaps then become 

N N N 

(vi/«({A})ivi/(2)(M)) =nrfi(A^.) i(oinc'i(A.)n^i(/^')io)^ 

j=i fc=i 1=1 

N 

= X{d,{^.M'^^'\{\}W\{^.})). (20) 

i=i 

The scalar product (\E'*^^''({A})|\E'*^^^({/i})) is of a special type namely, the state 
(\E''^^)({A}) I is an eigenstate while is not. This type of scalar product can 

be represented efficiently as a determinant [2S1 ES]- To obtain the normalized quench 
overlap Fx as given in ([7]) one should divide (!20|) by the norms of the states [SOj [311 [32] 
before and after the quench. The quench overlap is now defined up to a momentum 
dependent phase, which is related to the representation of the i?-matrix. 

So far we considered periodic boundary conditions, however we can easily generalize 
this by considering twisted boundary conditions. In this case the scalar product is given 
in [33]. Since the twist of the state after the quench does not appear explicitly, one has 
the freedom to chose the twists before and after the quench independently. 

2.2. Geometric quench overlap: Open boundary conditions 

The geometric quench setup has a natural generalization to the perhaps more physical 
case of open boundary conditions. 
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For an integrable system with open boundary conditions one should define the 
boundary conditions in terms of matrices /^_(A) and i^+(A). For XXZ-hke models the 
diagonal solution of the boundary Yang-Baxter equation [34j is given by 

with K_ and Kj^ describing the left and right boundaries respectively. Each is 

parametrized by ^-,^+. Following [35j we introduce 'double-row' monodromy matrices 
U- and W+: 

«_(A) = T(A)A-.(A)f (A) = ( ^-j;;) ) (22) 

U'j;(\) = r(A)'°At>(A)f (A)'" = ( ■^+<^^j ^+1^1 j (23) 

where 

f (A) = {-l)^alT'\-\)al (24) 

and to denotes matrix transposition in the auxiliary space. The operators 
•A~±{\),B±{\),C±{\) and V^i^X) are called the boundary operators in contrast with 
74(A), 5(A), C(A) and D{X) which are called the bulk operators. Eigenstates can be 
constructed using boundary operators 

N N 

(*±({A})| = (0|nC±(A,), |^±({A})) = X[B^{\m- (25) 

k=l k=l 

The action of boundary operators on the vacuum is similar to that of the bulk operators, 
expect for different vacuum eigenvalues. 

We now consider a geometric quench where we keep the left boundary fixed and 
move the right boundary from Li to L2. This allows us to write everything in terms 
of B-{\) and C_(A) operators. Equivalently we could fix the right boundary and write 
everything in terms of i3+(A) and C+(A). Similar to the case of periodic boundary 
conditions we factorize the monodromy matrix as 

W_(A) = r2(A)Ti(A)ir_(A)fi(A)f2(A) 

= T2{X)U.^i%{\) (26) 

where U-^i is the double-row monodromy of the smaller system with the same left 
boundary. We can now decompose the i3-(A) operator as a combination of bulk and 
boundary operators: 

B_{\) = (A2(A)i3_,i(A) + i?2(A)P_,i(A)) A2(-A) 

- (A2(A)^_,i(A) + 52(A)C_,i(A)) S2(A). (27) 

We write the initial state as 

N \ 

L(0|J]C_,i(Afe) 2(01, (28) 

k=l J 
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and we write the final state as 

JJ-B_,i(^fc)|0)i J ® I JjA2(^fc)A2(-//fc)|0)2 j + terms without overlap (29) 

^A:=l / \k=l / 

The terms without overlap are terms with no B-^i operator acting on |0)i and will 
not contribute when taking the scalar product with (!28l) . This results in the geometric 
quench overlap: 

N N N 

(vi/«({A})ivi/L^)(M))=n«2(^>2(-/i,)i(oinc4A.)n^-(/^oi^ 

i=l k=l 1=1 

M 

= U^2{^i>^{-f^,){^^-\{xm^ym))- (30) 

3=1 

For the XXZ spin chain, the determinant representation for the scalar product 
(^L^^({A})|^L^^({/i})) is computed in [361 |37] and is related to the overlap Fi of dH) 
by a normalization prefactor |36l[37j. Since ( l30l) does not depend on the right boundary 
condition, we can not only expand the system but also change the right boundary 
condition after the quench. 



3. Conclusions and perspectives 

We have shown how a sudden, global change in the geometrical support of an integrable 
model could be completely described within Bethe Ansatz, by providing an exact 
decomposition of the original (arbitrarily chosen) eigenstate into the basis of eigenstates 
after such a geometric quench. The idea is quite simple, but potentially very powerful: 
it applies directly to any model for which the Slavnov formula for wavefunction 
overlaps is explicitly known, and allows at least in principle to reconstruct the exact 
subsequent time evolution of wavefunctions, operator expectation values and time- 
dependent correlation functions via techniques similar to the ones used in the equilibrium 
case. On the practical side, since the energy of the states populated by the geometric 
quench can be obtained nonperturbatively, results can be accurately given for arbitrarily 
long times after the quench. Moreover, since the space of eigenstates is mapped out in 
detail via the Bethe Ansatz, extremely efficient truncations of the required Hilbert space 
sums can be implemented. We will explore a number of specific cases in forthcoming 
publications. 
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